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d ' Abstract 

We solve, by separation of variables, the problem of three anyons with a harmonic oscillator 
potential. The anyonic symmetry conditions from cyclic permutations are separable in our 
coordinates. The conditions from two-particle transpositions are not separable, but can be 
expressed as reflection symmetry conditions on the wave function and its normal derivative on 
the boundary of a circle. Thus the problem becomes one-dimensional. We solve this problem 
numerically by discretization. TV-point discretization with very small A'^ is often a good first 
approximation, on the other hand convergence as ^ cxd is sometimes very slow. 



1 Introduction 



Anyons are two-dimensional point particles with fractional statistics |ll|-|5|, which seem to be 
useful as models for the quasiparticle excitations in the fractional quantum Hall systems. For 
reviews see e.g. |^-10|. 



We use here what is known as the "anyonic" gauge, in which the vector potential associated 
with particle statistics vanishes, and the wave functions are many-valued, as follows. Go 
along a continuous curve in the A^-anyon configuration space, starting and ending with the 
same configuration A, but such that two particles are interchanged in the counterclockwise 
direction, while no other particles are encircled or interchanged. Then the value of a wave 
function ip must change continuously from t(^{A) to e'^ ip{A). The phase angle is a constant, 
independent of the configuration A and the wave function ip. It defines the statistics of the 
particles, it is for bosons, vr for fermions, and may take any real value for anyons. 

The energy spectrum for two anyons in a harmonic oscillator potential is easily found, 
and it interpolates continuously between the boson and fermion spectra when 9 varies Q. 
The energies vary linearly with 6, sometimes with a discontinuous change in slope at the 
Bose points. A more general class of quadratic Hamiltonians, including the case of a constant 
magnetic field, can be treated just as easily. 

The same problem for more than two anyons is much more difficult. Wu found a class 
of exact solutions that generalize the two-anyon solutions |lll| . More general exact solutions 



have been found, but all have energies that depend linearly on 9 [12-31, |]. Thus, even the 
three- anyon ground state close to Fermi statistics, which has an energy varying non- linearly 
with 9, is still not known analytically. Even though the exact ^-dependence is unknown for 
many energy levels, the spectral flow with 9 from the boson to the fermion spectrum is known 

Sen found an interesting "supersymmetry" in the three-anyon spectrum, excluding some of 



the states with linearly varying energy 23]. His supersymmetry transformation transforms 
bosons into fermions and vice versa, and more generally transforms 9 into tt — 9. 

The lowest part of the complete energy spectrum of three or four anyons in a harmonic 
oscillator potential has been calculated numerically |]32|-|3^. In the three-anyon problem, 
good analytical approximations to the wave functions corresponding to non-linear variation 



of energy are known [37|. More systematic approximation methods are perturbation theory. 



starting from the known boson and fermion spectra [38-40|, and the Hartree-Fock approxi 
mation |41 



One reason for the interest in the spectra of two- and three-particle systems in an external 
harmonic oscillator potential is that they can be used for calculating the second and third 
virial coefficients in the equation of state for an ideal gas of free particles, i.e. when there is 
no interaction apart from the statistics interaction. We will not discuss this problem here, 
but refer to the literature ||-|4|Jl|,|2|,||,||,|6| . 



We show here how to solve the three-anyon problem with a harmonic oscillator potential, 
by separation of variables in a suitable set of coordinates. The same method applies to a class 
of similar problems, as mentioned above. An anticlockwise cyclic interchange of the three 
particles gives a phase factor of e^^^ in the wave function, and this condition is compatible with 
the separation of variables. However, the interchange of only two particles gives a condition 
which can in general only be satisfied by a superposition of separated wave functions. Thus 
the separation is incomplete. We find general solutions numerically, but we hope that this 
method may also lead to some progress in the search for analytical solutions. 
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2 Formulation of the Problem 



It is a well known procedure to describe a particle in the plane by a complex coordinate z and 
its complex conjugate z, which may be scaled so as to become dimensionless. In the three- 
anyon problem the key to the separation of variables is a transformation from the particle 
coordinates 21,2:2,2:3 to the centre of mass coordinate Z and the relative coordinates u, f, 
defined by 

Z = (21 + 22 + 23) , 

u = (zi + rjZ2 + r/^23) , (1) 



--j= (^21 + r?22 + r/ 23 
1 

7i 



21 + ?7^22 + r/23 



Here r] = exp(2i7r/3) = ( — 1 + i\/3)/2 is a cube root of unity, with rf = ^ = rj and 
1 + r] + tf' = 0. Hence, 

Ti = (77(22 - 21) + 7/^(23 - 21 

V = ^[if{z2-zi) + i^{z^-zi)). (2) 

This coordinate transformation is a discrete Fourier transformation, which means that a 
cyclic interchange of particle positions, 

(^1, 2^2, 23) ^ Z2, 23) = (22, 23, 21) , (3) 

becomes diagonal, 

(Z, M, t;) {Z, u, v) = {Z, ifu, rjv) . (4) 

The interchange of particles 2 and 3 is just an interchange of u and v. Note that these two 
permutations generate the whole symmetric group S3. A similar treatment of permutation 
symmetry is known in nuclear physics p9| , [50| ] (we thank Alex Lande for pointing this out). 

Three particles in the plane define a triangle. The ratio s = (23 — 21)/ (22 — 21) is real when 
the triangle is degenerate so that the particles lie on a straight line. We define the orientation 
of a non-degenerate triangle as positive or negative depending on whether the imaginary part 
of s is positive or negative. We have that 

\u\ _ \r] + r]'^s\ _ Is + r/^l 
\v\ Ir/^ + rjsl \s + r]\ 

Hence |n| = \v\ if s is real, \u\ < \v\ if the orientation of the triangle is positive, and |n| > \v\ 
if the orientation is negative. 

The quantization of the centre of mass motion is trivial, and the problem we want to 
discuss here is the simultaneous diagonalization of the dimensionless relative Hamiltonian 
i^rei and relative angular momentum Ljei, defined by 

-f^rcl = - a a- - + UU + VV , 

ou ou ov ov 

d d _ d _ d 



Lrel = U^ + V- U — -V — 

ou ov ou ov 
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3 Separation of Variables 

The three-particle configuration is completely described by a total scale factor r > 0, a relative 
scale factor q >0, and two angles ipi and ip2 such that 



(6) 



These are the hyper spherical coordinates of Kilpatrick and Larsen [^], except that they use 
ifi ± tp2 instead of ^pi and ip2- We now have that 

^ 1 d d l+q^ ( l + q^ d ^ I ^ ^ 

2r^ dr dr 2r^ \ 9 dq dq q^ dtpi^ dip2'^ J 2 ' 

and 

Assume that the wave function of the relative motion is 

V' = V(r,g,^i,^2) = /(r)<7(9)e'(-''^^+''^^) . (9) 
Then the eigenvalue equation H^ei ip = Eip separates into an angular eigenvalue equation, 



(l + q^)[-^±q± + ^:l + k^]g = Xg, (10) 
y q dq aq q'^ 

with A as eigenvalue, and a radial equation. 



- I f"{r) - |- /'(r) + A + ^ j f^r) = Ef{r) . (11) 

A general wave function can be written as a linear combination of such separated wave 
functions. As shown in the next section, we need linear combinations, where A and j + k are 
constant but j — k varies, in order to satisfy the anyonic boundary conditions. 
In the usual way we find the solutions of the radial equation to be 

/(r)=r'^e-^'/2^a^r2™, (12) 

m=0 

with Ur = 0, 1, 2, . . ., with fi > —2 for normalizability, A = + 2), and 

E = 2 + f^ + 2nr . (13) 

The coefficients in the polynomial satisfy the recursion formula 

2(m + l)(m + 2 + ^) Um+i = {2m + 2 + fi - E) am ■ (14) 

The angular equation has two asymptotic solutions q^^ in the limit q ^ 0. We exclude 
the singular solution (for j = the singularity is logarithmic). In fact there is no reason for 
any singularity at q = 0, where the configuration is an equilateral triangle. This leaves a 
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solution which is unique up to a multiphcative constant, and with a convenient normahzation 
it may be written as 

g{q)=q\^\{l + q^rF{a,b;c--q^). (15) 
The constant k here may be chosen in one of two ways, 

K = | + l or «^ = -f' (16) 
giving two different representations of the same solution. The constants 

aJA^^., .= ^ + «. .= l + lil. (17) 

define the hypergeometric series 



m=0 



(c)m m\ 



where, e.g., (a)o = 1, (a)„+i = (a)„ (a + n). The convergence radius for this series is 1. There 
exist however other representations, 

g{q) = q\i\{l+q^Y-F{a,c-h-c-,^ 




= q\^\{l+q^Y-^Fyc-aMc;r^'j , (19) 

where the series converge for all real q. 

A useful parameter, as will be seen below, is the logarithmic derivative of 5 at g = 1, 

g{l) cF{a,b;c;-l) 

4 The Anyonic Boundary Conditions 

For three identical particles there is a six-fold identification of points in the Euclidean relative 
space. We will restrict the wave functions to the region < q < 1, which corresponds to 
all the positively oriented triangles, but is still a three- fold covering of the true configuration 
space. The boundary conditions defining the particles to be anyons are of two types, since 
there are two classes of non-trivial permutations. The first class contains the three-particle 
cyclic permutations, which leave q invariant. The second class contains the two-particle 
interchanges, which transform q into 1/q, and so give boundary conditions at q = 1. 

Consider first a continuous, counterclockwise and cyclic deformation of the configuration, 
as defined in eq. (^) , with no extra overall rotation of the triangle. It gives a phase factor e^'^ 
in the wave function, where 9 = vtt is the statistics parameter. If we deform while keeping all 
the time \u\ < \v\, the phase of v increases continuously from (p2 to ip2 + (27r/3), whereas the 
phase of u changes from ifi to (fi — (27r/3) + 2m'7r, where m' is any integer. The corresponding 
boundary condition on the wave function is, therefore, 

'4)[r,q,ipi- — + 2mTi,ip2 + —\=e 'il^{r,q,ipi,ip2) . (21) 
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j I h 2m'7r ] + k — = 2{n' + u)7r , (22) 



That is, 



for some integer n'. Since m' is an arbitrary integer, j must be an integer. Then 

k = j + 3{n + iy), (23) 

where n = n' — jm! is an arbitrary integer, and the eigenvalue of the relative angular momen- 
tum Lrei is 

(, = j ^k = 2j + 2,{n + v) . (24) 

These relations take care of the cyclic permutations of all three particles. What remains 
is only to take care of one of the three cases where two particles are interchanged, for example 
Z2 <-> 2:3, or equivalently, u ^ v. This is the same as q <-> 1/q and <-> ip2, if we define 
angles so that u = v corresponds to ipi = (p2- To be more precise, we consider a continuous 
interchange, with g = 1 at the beginning and end, and q < 1 during the interchange. The 
interchange should be anticlockwise, which means that we start with (pi > ip2, and end with 
(fi < (p2. There is one further restriction, that \(pi — (p2\ < (27r/3) when q = 1, meaning that 
the particle position zi must not be encircled. 

Thus, the boundary condition on ip at q = 1 is 

ip{r,l,ip2,ipi) = e''' ip{r,l,ipi,ip2) , (25) 

for < 1^1 — 1^2 < (27r/3). However, this condition is incomplete, because the general condition 
is 

'^l^ir, l/g,V32,¥'i) = e''^ i^{r,q,^i,y^2) ■ (26) 

Since the Schrodinger equation is second order in the q derivative, we need boundary con- 
ditions at q = 1 both for the wave function tp and its normal derivative ipq = dip/dq. The 
derivative condition is easily deduced, 

V'g(r, l,y?2,¥'i) = -e''^ i^g{r,l,^i,y^2) ■ (27) 

The boundary conditions for ip and 'i/'g can not be satisfied by a wave function which 
is separable in q, ipi and ip2- But we may quantize the relative angular momentum i, and 
according to eq. ( [2^ ) £ — = 2j + 3n is an integer, either even or odd. Let v' = u it n = 2m 
and u' = I' + 1 it n = 2m + 1, with m integer. Then 

j = Slrn+ — ] , k= - + 3{m + — ] . (28) 



Let gmiq) be the function g{q) as given by eq. ([l5|). Introducing ip = {ipi + ip2)/2, and 
summing over m, including an as yet undetermined coefficient Cm for each m, we get the 
following angular wave function, 

00 00 
niq,^i,^2)= E C'^5m(g)e^(^^'^^+'='^^)=e'^'^ C^5™(g) e"=^'(™+(^™^^-'^^) . (29) 



5 



It is natural to call U an anyonic spherical harmonic function, whenever it satisfies the anyonic 
boundary conditions. 

Define ^ = 3{ipi — <f2)- The two boundary conditions that must hold for < ^ < 27r are 

OD OO 

m=— OO m=— OO 

OO OO 

^ C^5;„(l)e'-« = -eH— Yl Cm9'm{l)e--^^ . (30) 

m=— OO m=— OO 

Recall that 5^(1) and g'^i^) depend on the three parameters fi = E — 2 — 2nr, £ and m. We 
claim that for each given i, the parameter which determines the energy E, may be adjusted 
so that the above boundary conditions have non-trivial solutions for the coefficients C^. For 
each £ there will be many solutions, possibly many with the same n, and this procedure should 
give the complete set of anyonic spherical harmonics. A complete proof would involve a study 
of the anyonic spherical harmonics, which we have not yet done, but wc can give a plausibility 
argument, which also indicates one possible way to solve the problem numerically. 

We may regard 7^ = Cmgmi^) and 7^ = Cmg'mi^) as the Fourier components of two 
functions 

00 

7(0= E C^-5m(l)e'"«, 

m=— 00 

00 

7(0= E Cr^g^e'"^, (31) 

m=— 00 

periodic in ^ with period 27r. We may restrict their regions of definition to be the interval 
[0,27r]. 

There is a natural scalar product between any two functions (p = ^(0 and x = x(0' with 
Fourier components <pm and Xm as above. 



icf>, x) = X(0 = 27r E (^rnT Xm ■ (32) 



r27r 
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Define the linear operator A by 

[A^](0 = e'(— '^'«^(27r-0, (33) 

for < ^ < 27r. Then A is Hermitean with respect to the natural scalar product, and A"^ = I, 
the identity operator. Note that ^ is a somewhat singular operator, unless v is an integer, 
since the factor e^i'^'^-'^ ^ extended by periodicity outside the interval [0, 27r], is discontinuous 
at every integer multiple of 27r. 

Define another operator B such that Bj = 7. It is also Hermitean, since it is diagonalized 
by the Fourier transformation, and its eigenvalues Pm = 5m(l)/5'm(l) are real. The two 
boundary conditions are simply 7 = Aj, Bj = —ABj, equivalent to two simultaneous 
eigenvalue equations, 

^7 = 7, {I + A)B{I + A)j = . (34) 

These two equations are compatible, since the operators A and B^ = (/ + A)B(I + A) 
commute, and B^ has a complete set of real eigenvalues, since it is Hermitean. If we select 



6 



one eigenvalue of 5+ on the subspace where A'j = 7, and set it equal to zero, this is one real 
equation for one real parameter /i. By equating different eigenvalues to zero we should obtain 
a complete set of solutions. 

There is one minor problem with this argument in that B is undefined if an eigenvalue 
Pm = g'mi^) / 9'm{^) is infinite. There is of course a similar argument using B^^ instead of B, 
but in principle it might happen that B and B~^ are undefined simultaneously. 



5 Numerical Solution 

To find a numerical solution for the coefficients Cm satisfying the boundary conditions in 
eq. (|30[), we must truncate to a finite number N of coefficients, and in order to use the fast 
Fourier transform, we choose A'^ as a power of 2. Then we impose eq. (^) at the N discrete 
points 

6 = ^^^^^, for A; = l,2,...,iV. (35) 

It is a remarkable feature of our method that for very small A^, say = 4, it can give 
many energy levels with good accuracy. This is so when the low Fourier components, i.e. with 
small m, dominate. On the other hand, the convergence as A^ ^ 00 is sometimes very slow. 
This is clearly related to the fact that the wave functions for non-integer u have non-integer 
power behaviour at ^ = 0, where two particles meet. Hence our approximations by finite 
Fourier series converge slowly. Note that eq. (l30|) at ^ = gives that 



(l_e--)7(0) = (l + e^'^-)7(0) = 0, (36) 

assuming that the wave function is continuous. Thus 7(0) = 0, except in the boson case, and 
7(0) = 0, except in the fermion case. For our numerical solutions these conditions at ^ = 
are not imposed and will hold only in the asymptotic limit A^ — > 00. 

To illustrate the convergence we have tabulated the approximate energy En as a function 
of A^, in three cases. Table 1 is for the state which becomes the bosonic ground state at 
V = Q. It is exactly known and has energy E = 2 + 'iv and angular momentum I = 'iv. 
Table 2 is for the state which becomes the fermionic ground state at u = \. It has angular 
momentum H. = —3 + 2>v and an energy depending on v in an unknown way. Table 3 is for 
the state connected to the fermionic ground state by the supersymmetry transformation of 



Sen |2^, 23 1 . It has angular momentum i = —2 + 'iv. The energy eigenvalues in Tables 2 and 
3 should be the same, except that v in one table corresponds to 1 — in the other. 

It appears from the tables that the leading correction term for finite A^ is of order N~'^'^ . 
Using two different A^ one may therefore extrapolate to A^ = cxd, and this improves the 
convergence considerably. Another point to note is that one may take advantage of the 
supersymmetry in order to get more accurate energy levels. 
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Table 1: The bosonic ground state energy 





u = 0.2 


1/ = 0.4 


= 0.6 


= 0.8 


-^exact 


2.6 


3.2 


3.8 


4.4 


N = 2 


2.3279 


2.972 08 


3.666 06815 


4.345 506 34 


N = 4 


2.3786 


3.06317 


3.743 586 44 


4.384 285 82 


N = 8 


2.4228 


3.11916 


3.775 664 77 


4.395 001 04 


N = 16 


2.4599 


3.152 74 


3.789423 68 


4.398 366 91 


N = S2 


2.4902 


3.172 56 


3.795 395 68 


4.399462 67 


N = 64 


2.5146 


3.18414 


3.797995 23 


4.399 822 87 


7 from fit to En = Ecxact + AN 


-7 




N = 2,4 


0.298 


0.7362 


1.247 39 


1.794 02 


N = 4,8 


0.321 


0.7592 


1.213 00 


1.652 37 


N = 8,16 


0.339 


0.7745 


1.202 21 


1.614 02 


N = 16, 32 


0.352 


0.7845 


1.199 78 


1.603 72 


N = 32, 64 


0.362 


0.7908 


1.199 56 


1.600 98 


£^oo from fit to En = 


Eoo + ^A^- 






AT = 2, 4 


2.5375 


3.186 09 


3.803 335 54 


4.403 375 54 


A^ = 4, 8 


2.5611 


3.194 70 


3.800 389 92 


4.400 275 74 


iV = 8, 16 


2.5759 


3.198 05 


3.800 028 70 


4.400 023 81 


TV = 16, 32 


2.5851 


3.199 31 


3.799 998 73 


4.400 002 07 


TV = 32, 64 


2.5909 


3.199 76 


3.799 998 91 


4.400 00018 
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Table 2: The fermionic ground state energy 





u = 0.2 


= 0.4 


u = 0.6 


ly = 0.8 


N 


= 2 


4.6455 


4.380 322 


4.181 3475 


4.048 5533 


N 


= 4 


4.6356 


4.379 612 


4.1818967 


4.047 9552 


N 


= 8 


4.6446 


4.394153 


4.189 5454 


4.049 7399 


N 


= 16 


4.6552 


4.405 332 


4.194 0494 


4.050 6380 


N 


= 32 


4.6646 


4.412 454 


4.196 2485 


4.050 9943 


N 


= 64 


4.6723 


4.416 732 


4.197 2586 


4.051 1243 


N 


= 128 


4.6786 


4.419 241 


4.197 7105 


4.051 1700 


7 from fit to £^iv 


= Eoo + AN-^ 






N 


= 8, 16, 32 


0.171 


0.650 


1.034 


1.334 


N 


= 16, 32, 64 


0.276 


0.735 


1.122 


1.454 


N 


= 32, 64, 128 


0.319 


0.770 


1.160 


1.512 


from fit to En = E^c + AN'^ 


N 


= 8, 16, 32 


4.7394 


4.424 958 


4.198 3467 


4.051 2286 


N 


= 16, 32, 64 


4.7092 


4.423 167 


4.198 1166 


4.051 1991 


N 


= 32, 64, 128 


4.7038 


4.422 799 


4.198 0764 


4.051 1946 


Eoo from fit to En = E^o + AN''^'' 


N 


= 2,4 


4.6046 


4.378 652 


4.182 3200 


4.0476609 


N 


= 4,8 


4.6726 


4.413 775 


4.195 4408 


4.050 6184 


N 


= 8, 16 


4.6883 


4.420 416 


4.197 5209 


4.051 0800 


N 


= 16, 32 


4.6940 


4.422 065 


4.197 9435 


4.051 1697 


N 


= 32, 64 


4.6966 


4.422 505 


4.198 0371 


4.051 1884 


N 


= 6i. 128 


1.(3981 


1. 122 582 


1.198 0588 


1.051 1924 
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Table 3: The supersymmetric partner of the fermionic ground state 





u = 0.2 


u = 0.4 


u = 0.6 


ly = 0.8 




= 2 


3.8436 


3.954 56 


4.249 4885 


4.627 147 52 


N 


= 4 


3.8814 


4.051 14 


4.348 6537 


4.678 013 75 


N 


= 8 


3.9144 


4.11177 


4.391 1581 


4.692 955 72 


N 


= 16 


3.9423 


4.147 84 


4.409 1615 


4.697 703 46 


N 


= 32 


3.9653 


4.168 98 


4.416 8475 


4.699 244 71 


N 


= 64 


3.9840 


4.18127 


4.420 1526 


4.699 74811 


7 from fit to En = Eoc + AN'^ 


N 


= 2, 4, 8 


0.194 


0.672 


1.222 


1.767 


N 


= 4, 8, 16 


0.246 


0.749 


1.239 


1.654 


N 


= 8, 16, 32 


0.275 


0.772 


1.228 


1.623 


N 


= 16, 32, 64 


0.299 


0.782 


1.218 


1.614 


-Boo from fit to E^ = E^^ 


+ AiV-T 






N 


= 2, 4, 8 


4.1441 


4.21400 


4.423 0430 


4.699 170 53 


N 


= 4, 8, 16 


4.0925 


4.200 86 


4.422 3905 


4.699 91463 


N 


= 8, 16, 32 


4.0749 


4.198 87 


4.422 5732 


4.699 985 54 


N 


= 16, 32, 64 


4.0652 


4.198 36 


4.422 6461 


4.699 992 28 


Eoo from fit to £^Ar = E^o 


+ AN-^" 






N 


= 2,4 


3.9997 


4.18146 


4.425 0877 


4.703 053 33 


N 


= 4, 8 


4.0178 


4.193 57 


4.423 9194 


4.700 311 10 


N 


= 8, 16 


4.0295 


4.196 52 


4.423 0381 


4.700 040 61 


N 


= 16, 32 


4.0374 


4.19749 


4.422 7717 


4.700 003 41 


N 


= 32, 61 


4.012G 


1.197 86 


4.422 7001 


4.699 995 92 
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